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ABSTRACT 

We consider the Raylcigh-Taylor instability in the early evolution of the rarefied 
radio bubbles (cavities) observed in many cooling-flow clusters of galaxies. The top of 
a bubble becomes prone to the Rayleigh-Taylor instability as the bubble rises through 
■^j- ' the intra-cluster medium (ICM). We show that while the jet is powering the inflation, 

C^) . the deceleration of the bubble-ICM interface is able to reverse the Rayleigh-Taylor 

JO ' instability criterion. In addition, the inflation introduces a drag effect which increases 

substantially the instability growth time. The combined action of these two effects 
considerably delays the onset of the instability. Later on, when the magnitude of the 
, deceleration drops or the jet fades, the Rayleigh- Taylor and the Kelvin Hclmholtz 

instabilities set in and eventually disrupt the bubble. We conclude that the initial 
. deceleration and drag, albeit unable to prevent the disruption of a bubble, may signif- 
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1 INTRODUCTION 

The recent high spatial-resolution observations show that many clusters and groups of galaxies harbour bubbles (cavities) 
devoid of X-ray emission. Examples include Perseus (Fabian et al. 2000), Hydra A (McNamara et al. 2000), Abell 2052 
(Blanton et al. 2001), Abell 2597 (McNamara et al. 2001), MKW 3s (Mazzotta et al. 2002), HCG 62 (Vrtilek et al. 2002), 
Abell 4059 (Heinz et al. 2002) Abell 478 (Sun et al. 2003), and MS 0735.6+7421 (McNamara et al. 2005); see also Birzan 
et al. (2004) for a systematic study. These bubbles are very low density regions inflated by the jets launched by the active 
galactic nuclei (AGN) sitting at the centres of these clusters. The lack of strong shocks at the rims of these bubbles shows that 
their inflation is quite gentle, and the bubbles are in approximate pressure equilibrium with the surrounding environment. On 
account of their lower density with respect to the outer environment, the bubbles rise buoyantly and sweep up some of the 
intra-cluster medium (ICM), as evident by the enhanced emission at their rims. 

It would seem that the top of these bubbles is prone to the Rayleigh-Taylor Instability (hereafter, RTI). Indeed, the 
gravitational acceleration of the cluster pushes downwards, and the dense swept-up gas shell lies on the top of the thinner 
one inside the bubble itself. This RTI would be able to tear the bubble apart in few characteristic e-folding times 
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where \g\ the modulus of the cluster gravitational acceleration and A/j is the wavelength of the RT perturbation. The char- 
acteristic time (1) is about one order of magnitude less than the estimated ages of most bubbles (Birzan et al. 2004). This 
observation has prompted many authors (e.g. Briiggen and Kaiser 2001; Kaiser et al. 2005; Reynolds et al. 2005; Jones and 
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De Young 2005) to invoke an ordered magnetic field at the edge of the bubble to stabilise the bubble against the RTI. This 
scenario, however, may lead to some difficulties: if the magnetic field is ordered on a scale of the same order as the size of 
the bubble, its effect should be highly non-isotropic: it should be able to suppress very efficiently the RTI modes along the 
magnetic field lines, but not perpendicular to them. This is because the magnetic field cannot exert stress countering the 
instability a direction perpendicular to the field lines. If the field is tangled on scales much smaller than the bubble size, then 
it will have no effect on the large-wavelength modes which are those that tear the bubble apart. 

In this paper we argue that a magnetic field is actually not necessary to explain the long lives of the radio bubbles. 
Extending the work by Soker et al. (2002) we shall show that during the first stages of its inflation a bubble is stable. In the 
early phase the inflation of the bubble is powered by a jet launched by the central AGN. While the jet is active the bubble 
radius grows with time as Rt(t) oc t 3/l5 , implying that the acceleration Rb is negative. On the top of the bubble, then, this 
deceleration overcomes gravity, and the configuration of the bubble is RT-stable. Besides, the inflation introduces a drag 
effect which greatly reduces the instabilities growth rate. Only later, when \Rb\ decreases and the gravitational acceleration 
prevails, or the jet fades, does the instability set in, tearing the bubble apart. We show that the delay of the RTI onset is 
more than enough to explain the survival of old bubbles. 

This paper is laid out as follows; in § 2 we discuss the properties of the bubbles, and sketch the physics developed in 
this paper. In § 3 we derive the equations of the RTI in a bubble actively inflated by a jet. In § 4 these equations are solved 
numerically, and the results are discussed in § 5. The last stage of the life of a bubble, after the extinction of the jet, is 
addressed in § 6. We summarise in § 7. Readers not interested in the mathematical details may harmlessly skip § 3 and § 4, 
jumping directly from § 2 to § 5. 



2 THE RAYLEIGH— TAYLOR INSTABILITY IN BUBBLES 

In this Section we refine the treatment by Soker et al. (2002) of the RTI of a bubble inflated by an AGN jet. The jet stops 
at a distance r c from the cluster core, and hence it becomes the expansion centre of the bubble. The injected jet's material 
passes through a strong shock and forms a very hot bubble, which later on keeps in approximate pressure equilibrium with 
its surroundings. On account of its low density, the bubble rises buoyantly into the ICM; during this process the expanding 
bubble sweeps up the ICM to form a thick shell of denser gas ahead of itself, as evident from the local enhanced X-ray 
luminosity. The X-ray observations also show a lack of strong shocks ahead of the bubble, so the swept-up shell and the gas 
inside the bubble move about with the same velocity. Because the shell density p s is much higher than the density pb inside 
the bubble, and since the gravitational acceleration is directed downward, the flow at the bubble's top would be RT unstable 
if only gravity is considered. This simple picture, however, overlooks the bubble's inflation process. We therefore consider a 
local frame of reference attached to a point P on the surface of the expanding bubble, and moving with it. A test particle in 
the non-inertial frame attached to P feels the radial (i.e., normal to the bubble's surface) acceleration 

g = -\g\ cos8-R b (t), (2) 

where g is the cluster's gravitational acceleration, the angle between the orthogonal to the bubble surface in P and the 
radial direction of the cluster, and R b (t) is the bubble radius at the instant t. If the jet supplies energy at the constant rate 
of E, then 

R b (t) = a (Et 3 / Pa y /J (3) 

(Castor et al. 1975; Dokuchaev 2002), where p a is the density of the ambient ICM. The dimensionless parameter a depends 
on the equation of state (EOS) of the gas inside the bubble: for a non-relativistic EOS a = 0.929, while for a relativistic EOS 
a = 0.793 (Dokuchaev 2002). For bubbles in cooling flow clusters the EOS is most likely to be somewhere in between these 
two extremes. It is worth noticing that Equation (3) assumes a uniform ambient density, which is not true for the stratified 
atmosphere of a cluster the bubble is embedded in. Nevertheless, in the following we shall adopt Equation (3) since it is but 
a secondary source of error. Moreover, if the cluster's negative density gradient is considered, then the bubble's top has the 
largest acceleration, making it more stable. Since the bubble's top front is the least stable region, considering the density 
gradient will make the stabilising effect studied here more favourable. 

The normal acceleration of a test particle residing on the bubble-ICM boundary in the non-inertial frame is therefore 

g = -\g\ cos6+±^fi. (4) 

Near the bottom of the bubble 9 ~ tt, so g is always positive, meaning that the cluster's gravitational acceleration and R b 
have the same direction, and the total acceleration pulls towards the cluster centre. In this segment, the bubble is RT stable, 
because the lighter plasma (inside the bubble) lies above the denser one. Near the top of the bubble (9 ~ 0) the acceleration (4) 
may be either positive or negative. At early stages (t small), Rb is large enough to overcome gravity, thus making g positive. 
The acceleration pulls upwards, and on account of the density stratification even this segment is RT-stable, without the need 
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of any magnetic field (as pointed out by Soker et al. 2002). On top of this, as discussed in § 3, the inflation introduces a 
drag effect on the perturbation, akin to the Hubble drag affecting the perturbations growth in an expanding Universe. Both 
in the cosmological context and here, this drag reduces the perturbations growth rate, further delaying the RTI onset. Later 
on, however the gravitational pull will prevail, and the RTI commence. The delay introduced by the early contribution of R b 
and by the drag term is more than enough to lengthen the life time of the bubble to a significant amount. The mathematical 
details justifying this statement are worked out in the following Section. Readers not particularly eager to afford them may 
skip § 3 and § 4 and go directly to Section 5. 



3 RAYLEIGH-TAYLOR INSTABILITY OF AN EXPANDING BUBBLE 

In this Section we derive the equations for the Rayleigh-Taylor instability of an expanding bubble. As customary in any 
perturbation analysis, it is convenient to decouple the "background" motion of the bubble (i.e., its expansion and uprising) 
from the motion associated to the true instability. To this aim we introduce a frame of reference moving with the expanding 
bubble. We consider an inertial frame of reference K, the origin of which coincides with the cluster's centre. In this frame of 
reference the bubble's centre O has coordinate r b (t). We define a second frame K" with its axes parallel to those of K, and 
centred on r b . Finally, we define the non-inertial frame K' with the same centre as K" , which co-moves with the expanding 
bubble: in other words, the coordinates of any point expanding with the bubble are constant in the frame K' . The relation 
between the coordinates x" and r', referred to the frames K" and K' respectively, is 

x" — a(t) r , (5) 

where (from Equation (3)) we define the expansion scale factor 

a(t) = R b (t)/R b (t ) = (t/t ) 3/5 , (6) 

where to is an arbitrary reference time. Putting together all the previous equations, the relation between the inertial frame 
K and the co-moving frame K' is 

x = r b (t) +a(t) r. (7) 

This formula allows us to write down a set of hydrodynamic equations in the co-moving frame K' . The advantage of this 
approach is that the unperturbed equations in K' are time-independent, allowing a more direct approach to the subsequent 
perturbation analysis. 

We write the usual gas dynamic equations in the inertial frame K (e.g. Landau and Lifshits 1987) 

^+V-(p«)=0 (8a) 

P^ + VP-pg = 0. (8b) 

From the coordinate transformation (7) and the equality t = t' (K and K' have the same time units), we retrieve the relations 
between the derivatives in K and K': 

d d Vb + ar' , , . 

di = dF-^^- V (9a) 

V = i V'. (9b) 
a 

In the first equality we assume that during the relevant evolutionary phases for our study the buoyancy velocity v b = r b is 
negligible in comparison to the expansion velocity a r' . Plugging these into Equations (8), after some algebra, we get the 
following set of equations in the co-moving frame K'\ 

^+V-V«') = (10a) 

''(^ + r'+ 2 H + U v ' PW9 H (10b) 

where we have introduced the new co-moving variables 

p = p a 3 (U a ) 
' v-ax' 

v — (lib) 

a 

P' = Pa 2 , (11c) 
and d! /dt' = d t i +v' • V' is the convective derivative in the co-moving frame. 
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3.1 The Unperturbed Configuration 

As noted in the discussion following Equation (4), the most prone segment to the RTt is the top of the bubble. Where 
necessary, therefore, we shall limit our analysis to this portion of the bubble. Before studying the perturbations, it is necessary 
to define a "background" unperturbed configuration. In the present case it is, obviously, the configuration of an expanding 
bubble. This is described by the set of equations (10) with v = 0. The co-moving continuity equation (10a) yields d t > p = 0, 
i.e. p is constant with time; the density p in the inertial frame is p oc a -3 , i.e. the mass contained within any given volume 
expanding with the bubble is constant, as expected. 

As our X-ray observations show, there is no signature of strong shocks between the bubble interior and the overlying 
swept-up shell. For this reason, we assume that the shell and the adjacent segment of the bubble interior move with the 
same speed, and the only difference between them is in density. We denote with R the co-moving coordinate of the interface 
between the thick shell and the bubble interior. The unperturbed density profile undergoes a sharp discontinuity across the 
edge of the bubble, passing from the low value pt inside the bubble (r 1 < R) to the ambient value p a outside it (r' > R). We 
put p s ~ p a , i.e. we neglect the density enhancement ahead of the bubble due to the swept-up material, since its effect is 
small. Our unperturbed co-moving density profile is then 



pa r' > R 

pt r' < R. 



(12) 



The co-moving unperturbed static Euler equation 

VP' =p'(g-a r') , (13) 

can be projected in the radial and tangent direction with respect to the co-moving frame K': 

f)P' 

a — = p> {gr ,-ar\ (14) 

V h P' = p g h , (15) 

where g r i and gh are respectively the projections of the cluster gravitational acceleration parallel and perpendicular to the 
direction of the bubble's radius. It is quite apparent that the term gh is non-zero if the point P is not aligned with the line 
joining the bubble and the cluster centres. On the top of the bubble, which concerns us most, gh ~ 0, so VhP ~ 0, and the 
pressure gradient is purely radial. Here, Equation (14) defines a local hydrostatic equilibrium in the effective gravity 

g r , -tir' = - [\g\ + a(t) r'] = -g. (16) 

The overall effective gravity is smaller than expected from the cluster's contribution alone, since a(t) < 0. 



3.2 Lagrangian Perturbation Analysis 

Having determined the properties of the equilibrium solution, we pass to study the perturbation to the inflating bubble. For 
the sake of simplicity, from now on we omit all the primes referred to the co-moving frame K' . It should be clear, however, 
that all the variables hereafter refer to that frame. 

The perturbation analysis is most easily performed in the Lagrangian formalism (see e.g. Shapiro and Teukolsky 1983 
for a nice introduction). Let £(x,t) be the displacement of a fluid element due to the perturbation. The Lagrangian density 
perturbation Ap, i.e. the perturbation measured in a frame of reference moving with the unperturbed flow, is 

Ap=-pV-£. (17) 

In the context of the study of the RTI, it is customary to assume that the perturbation is incompressible. This assumption 
holds if the RTI growth time (Equation (1)) is shorter than the time 

*— ~ A/c * ~ L0 x 106 yr (-fe) (icTTev) • (18) 

taken by a sound wave to cross the perturbation itself. The rationale of this request is that if iaound is very short, then the 
sound waves are very effective in keeping the pressure equilibrium within the portions of the perturbation. In our case the 
incompressibility requirement t-RT <C tsound does not seem to hold. However, the compressibility affects little the growth rate 
of the RT perturbation when the density ratio pa/pt is large, as occurs in the present case (see e.g Baker 1983; Livescu 2004). 
Therefore, it is safe to adopt the incompressibility assumption Ap = 0: from Equation (17), this is equivalent to 

V-£ = 0. (19) 

The a-th component of the perturbed Euler equation reads 

d 2 £ a . „ a d£ a 



+ 2 - 

dt 2 a dt 



+ V a (AP)- (V S P ) Slats =0, (20) 
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where the pressure gradient may be expressed in terms of the unperturbed Euler equation VP = p g. It is convenient to 
decompose the displacement vector £ to a radial component £ and a tangential (or horizontal) component £h'- 



The condition of incompressibility (19) then reads 



— — ( 2 A 

r 2 dr A r 



£ 



V h -$ h = 0. 



(21) 



(22) 



(23) 



If we may neglect the radius of curvature of the bubble, this simplifies to 

|^ + Vh • £ h ~ 0. 

or 

This approximation is permitted so long as the size A of a perturbation (which is stretched by the bubble inflation itself) 
is much smaller than bubble radius Rt- Although we are also interested in perturbation with A ~ Rb, we will use this 
approximation throughout the rest of the paper, as it is adequate enough for our goal of showing that no stabilising magnetic 
fields are needed to explain the survival of old bubbles. 

We take advantage of the fact that the unperturbed quantities only depend on the radial coordinate r to separate out 
the dependence on the horizontal coordinates: 

£ = £(r,t) e ikh ' x (24a) 

£h = £h{r,t)B ih ^ x (24b) 

AP = AP{r,t) e lkh X , (24c) 

where kh is the horizontal (i.e., tangential) co-moving wavenumber vector. The perturbed Euler equation (20) reads, compo- 
nent by component, 

d 2 



dt 2 



£ 
£„ 



d 


' £ ' 


l 

+ - 


dAP/dr 


+ E1 

a 


d£/dr 


dt 




a 


i k h AP 


i k h £ 



= 0, 



where g is defined by Equation (16). We also have the (approximated) condition of incompressibility 

d£ 



dr 



+ i k h ■ £ h = 0. 



(25) 



(26) 



Dotting the horizontal component of the Euler equation by kh and eliminating from there £h ■ kh via the incompressibility 
condition, we retrieve 



AP = - P {g£ + ^ 



d*_d£ +2 & d_d£ 
dt 2 dr a dt dr 



Inserting this into the radial component, after some algebra we find 

d 



l 2 
P K h 



d 2 £ +2 dd£ 
dt 2 a dt 



dr 



(f_d£ 2 hd_ d£ 
dt 2 dr a dt dr 



- ( Pa - p b ) 6(r -R)i^£ = 0, 
a 



(27) 



(28) 



where we have neglected the spatial variation of the effective gravity g on the scales we are concerned with; the delta term 
stems from the radial derivative of the unperturbed discontinuous density profile (12). 

We further simplify the RTI criterion given by Equation (28). In our frame of reference the background is static, so 
d/dt = d/dt. Equation (28) must be solved separately inside and outside the bubble, i.e. for r § R. The two branches are not 
independent of each other, but are linked by the jump condition imposed by the delta function in Equation (28) . For r 7^ R 
Equation (28) reduces to 



1.2 

" h 



d 2 t £ + 2-d t £ 
a 



d t t d rr £ + 2 - d t d rr £ 
a 



0. 



For this equation we look for solutions in the separate form 

£(r,t) = cj>(t) C(r) 

We find 



dr 2 kh C 



0. 



so Equation (29) has the piecewise solution 



£{r,t) 



4>(t) e- kh (r - R) 
4>(t) e fe " (r ~ R) 



r>R 
r <R. 



(29) 



(30) 



(31) 



(32) 



The function <j>{t) is determined by the jump condition across the surface of discontinuity r = R. We integrate Equation (28) 
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over the thin layer [R — e, R + e], where e — > 0. After some algebra we derive 



lim 



0(t)(p Q - P6 ), (33) 

K-e 



where on the right-hand side we have used the identity </>(t) = £(r = R,t). 

We now may insert the solution (32) into the jump condition (33) to get our equation for the evolution of the Lagrangian 
perturbation: 

d 2 4> a dtp gk h p a - p b 

-JTT+2--J- ; 4> = 0, (34) 

dt z a dt a p a + pb 

where the effective acceleration g must be evaluated on the top of the bubble. 

Note that if a = constant this formula would reduce to the familiar RTI equation 

~Tl2 = \9\ k h — — 9, (35) 

dV pa + Pb 

the growing solution of which is the exponential cf> oc e n with 

n 2 = Ifll k h E?lZ£±, (36) 

Pa + Pb 

If p a > Pb, i.e. the denser fluid lies on top of the lighter one, then n 2 > 0: the perturbations grows exponentially as in the 
standard scenario of the RTI. 

In the case of the inflating bubble, it is convenient to express the coefficients of Equation (34) terms of physical (rather 
than co-moving) variables. The effective gravity on the top of the bubble is 

g = \g(r b + R b )\+R b (t), (37) 

The wavelength Xh = 2 Tr/ku in all our equations is a co-moving quantity, so its magnitude does not have a precise meaning per 
se, but only when it forms a dimensionless quantity with the curvature radius Rb of the bubble. The dimensionless quantity 

x = k h {t)R b (t), (38) 

i.e. the ratio between the bubble circumference and the perturbation's physical wavelength is frozen in time, since the pertur- 
bation is stretched by the expansion of the bubble. In the last equation we have written fc^(t) to remark that here this wave 
number refers to a wave in physical (i.e. non-comoving) coordinates. In order to be consistent with our planar approximation, 
the factor h must be large. Finally, since the bubble's density is always much lower than the ambient density 

Pb<.pa, (39) 

we write down the final expression of the perturbation equation in the form 

<F±,e}_d±_ \g\+R\{ t) 
dt 2 5 t dt * R b (t) 

where we used Equation (6) to substitute for a/a. It is worth to remember that <j> is a co-moving variable: its link with the 
"physical" amplitude <^ phys (£), i.e. the amplitude expressed in non-comoving coordinates is given by the relation 

<i>{t) = <f> phys (t)/R b (t). (41) 

The qualitative properties of Equation (40) are briefly discussed in the next Subsection. 



+i7T- ^ yi ,;:r ^o. (40) 



3.3 Properties of the Rayleigh— Taylor Instability of an Inflating Bubble 

Equation (40) differs from the standard equation for the Rayleigh-Taylor instability in three respects. 

(i) The term proportional to dcj>/dt. This is essentially due to the fact that the unperturbed bubble is inflating; this term is 
therefore akin to the "Hubble drag" in the equations for the evolution of cosmological perturbations in an expanding Universe 
(see e.g. Peacock 1999, page 470). 

(ii) The total acceleration g is made up by two contributions, the gravitational acceleration and the inflation deceleration. 
We have already illustrated how this may affect the onset of the RTI (see the discussion of Equation (4) in § 2). 

(iii) The last term in Equation (40) is proportional to -R^ 1 . Again, this is an effect of the co-moving frame: the factor R^ 1 
stems from the fact that on account of the stretching due to the inflation all the co-moving lengths are larger by a factor 
oc R b (t) than their non-comoving counterparts. 

It is also instructive to give some insight in the qualitative behaviour of the solution of Equation (40). The coefficients 
of Equation (40) are singular for t — > 0, which is expected because the inflation starts at t = 0. Since R b < 0, at early stages 
the acceleration is dominated by inflation, and Equation (40) can be approximated by 

I^ + l.i (42) 

dt 2 + 5 t dt + ^ * ¥ °- (42) 
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This equation admits analytical solution of the form <j> oc t v , where the complex index 

V = -w±JoVl-2A> i . (43) 

is solution to the quadratic equation 

V 2 + \V + A x = 0. (44) 

The wavenumber of a perturbation is clearly limited by the size of the bubble: since yc 1, the argument of the square root 
is always negative: rearranging the real and imaginary part of f, 

<(>(t) = {t/t )- 1/10 {Acos[uj log(t/t )] + Bsin[o; log(t/*o)], } (45) 
where to is an arbitrary time and 

u> = & |24^- 1| 1/2 . (46) 

Apart from the oscillating factor, the co-moving perturbation amplitude <$> decays as showing that the bubble is 

Rayleigh-Taylor stable. In non-comoving coordinates (see Equation (41)) this the amplitude grows as 4> p h_ ya oc t 1 ^ 2 , i.e. less 
than the bubble radius Rt oc t 3//5 . The acceleration is inflation-dominated, and is directed upwards: this and the relative 
stratification of the bubble and the overlying medium imply that the configuration is stable, as found. Note that the "Hubble 
drag" term in Equation (40) plays a fundamental role in this decay. Without this term the early-time equation (42) would 
read 

i\ 
dt 2 

Again, this admits solutions in the form 

<t>(t) = {t/t ) 1/2 {Acos[uj' log(V*o)] + Bsm[uj' log(t/*o)]} , (48) 

with 

w'= i |24x-25| 1/2 . (49) 

In this case (apart from a slowly oscillating term) the amplitude increases, albeit at the somewhat slow pace <j> oc t 1 ' 2 . 
This apparently strange result is correct, since the energy of the oscillations diminishes, as the following argument shows. 
Equation (47) is formally the equation of a harmonic oscillator with a time-dependent square angular frequency 

Vi 2 = A J. (50) 

It is well-known from analytical mechanics (see e.g. Arnold 1978) that the ratio 

I = E/Q. (51) 

between the energy E ~ Q 2 (j> 2 of the oscillator and its angular frequency Q is an adiabatic invariant. Therefore, since Q oc t _1 , 
then cf> oc t 1 / 2 and 

ExT 1 , (52) 

which shows that the energy of the oscillations decreases. 



^ + A>4~0. (47) 



4 NUMERICAL RESULTS 



In this Section we solve numerically Equation (40), which describes the evolution of the co-moving Lagrangian amplitude 
4> (Equation (41)) of a Rayleigh-Taylor perturbation in a non-inertial frame of reference moving with the expansion of the 
bubble, given by Equation (3). 

We refer our numerical calculations to the cluster A 2052: for the ambient density profile we adopt the f3 - model 

Pc 



p{r) 



(53) 



[1 + (r/r c ) 2 fP /2 

with p c = 5.7 x 10~ 26 g cm" 3 , r c = 36.4 kpc and /3 = 0.56, which provides an approximate fit to the deprojected density 
profile given by Blanton et al. (2001). 

For the gravitational acceleration we assume the Navarro et al. (1997) profile calculated by Zakamska and Narayan (2003) 
on the same data by Blanton et al. (2001): 

r/r s 



(r) = 2 g (r/r.Y 



log(l + r/r s ) 



1 + r/r s 



(54) 



where r s = 340 kpc and go = 2.8 x 10 cm s 
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The radius of the bubble (see also Equation (3) and Soker et al. 2002 ) is 

R ^- 7 - 8 kpC ( l0^ ergs-0 (lO^r) ( l^gcm- J ' (55) 
In our calculation we must also take into account the buoyancy of the bubble: the drag force limits the bubble rise velocity 
to the uniform terminal velocity 

v , = 8 \Ajm (56) 

where (following Churazov et al. 2001) we set the drag coefficient Cd ~ 0.75. This velocity is slow enough to allow us to 
simply update at every integration step the values of \g\ and p a with the values of the new environment in which the bubble 
is embedded. 

The choice of the initial integration time to an d of the initial condition are not fully trivial, on account of the singularity 
of the equation for t = 0. We choose an "early" value of to, meaning that at this time the inflation acceleration Rb is still 
dominant over gravity. At this time we fix the value of <j>, the ratio between the physical (i.e., non-co-moving) amplitude of 
the perturbation and the radius of the bubble (see Equation (41)): 

4>{t = t ) = 4>a. (57a) 

Differentiating Equation (41) we get 

0(t = to) = ^o = 0j£-fi); (57b) 

the derivative (j> measures the relative magnitude of the growth time scales of the perturbation t p = p hys/0phys and tb — 
Rb/Rb = 5to/3 of the bubble. We have integrated numerically Equation (40), with a fifth-order adaptive Runge-Kutta 
numerical algorithm (Press et al. 1992). We have explored its solutions for different values of the initial conditions (57). We 
have stopped (quite arbitrarily) the integration at t max = 10 s yr after the inflation starts. This time is few times larger than 
the life times of the radio bubbles estimated by Birzan et al. (2004). Of course, we have ascertained that our result do not 
depend on our choice of t max . 



5 DISCUSSION OF THE NUMERICAL RESULTS 

In Figure 1 we plot the evolution of Rb, Rb and Rb for the specific case of a bubble released at r c — 5 kpc from the cluster 
centre. In Figures 2, 3, 4, 5, 6 and 7 we plot the evolution of the co-moving amplitude cj>{t) for different parameters: 

• the initial perturbation amplitude (f>o (see Equation (57a)); 

• the ratio between the growth time t p of the perturbation and of the bubble tb (see Equation (57b); 

• the initial time to at which the integration starts 

• the distance r c from the cluster centre at which the bubble is initially inflated; 

• the factor x — kh Rb, i.e. the ratio between the circumference of the bubble and the horizontal wavelength of the physical 
RT perturbation. 

Finally, in Figure 8 we have plotted the modulus of the relative velocity of the RT perturbation with respect to the expanding 
front of the bubble. In all our calculations we have kept fixed the value of the energy injection rate E = 10 44 erg s _1 . 

In Figure (2) we show the effect of the different terms in Equation (40). The dashed line plots the modulus \<j)\ of the 
solution of Equation (40) in which the effects of the inflationary deceleration are neglected: 

w-"m>*-' y - (58) 

here inflation is only taken into account as a progressive growth of the bubble's radius Rb- The dot-dashed line plots the 
solution of Equation (40) in which the inflationary deceleration is taken into account, but where the "Hubble drag" term (i.e., 
the term proportional to <f>) is neglected: 

Finally, the solid line plots the solution of Equation (40) as it stands. The ICM density profile is given by Equation (53) and 
the gravity profile by Equation (54), and refer to the cluster A 2052. 

In all cases the amplitude grows, leading to the eventual bubble disruption by the RTI. At the "equivalence time" 
t eq ~ 2 x 10 7 yr when the moduli of the gravitational acceleration and the inflation deceleration are equal, the amplitude 



1 This formula assumes that the interior of the bubble is much thinner than the surrounding environment. 
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\4>\ is still small in the case described by Equation (40) (solid line), but is already well-developed in the case described by 
Equation (58) (dashed line) and in the case described by Equation (59) (dot-dashed line); in this latter case the perturbation 
oscillates, and on account of the adiabatic invariance of the ratio between the oscillations' energy and frequency, the amplitude 
increases as <j> oc t 1 ' 2 (see the discussion after Equation (47)). At t cq the magnitude of \Rb\ falls below the value of gravity, 
and the overall acceleration changes sign; the RTI sets in and after few t eq the bubble is torn apart (but see also the next 
Section for more important details). 

It may seem a little puzzling that the inclusion of the deceleration Rb may shorten the life time of the bubble. This 
behaviour depends on the magnitude of the quantity n 2 = >c \g\/Rb, as we show with the following argument. If we ignore 
the time dependence of n 2 , then the solution of (58) grows exponentially: 

<^oce nt . (60) 

On the other hand, Equation (59) at small times when inflation dominates has solution 

(j><xt 1/2 , (61) 

as discussed in § 3.3. In the extreme case n 2 1, then the solution (60) grows very fast, while the solution (61) is slower. 
At the opposite extreme n 2 < 1 the solution (60) keeps limited for t < n -1 , and is rapidly outgrown by the solution (61). 
Besides, from t eq on, the oscillations in Equation (59) become unstable. By this time they have developed to a fairly large 
amplitude, which promptly blows up as soon as the instability sets in. Since the factor n 2 = x\g\/Rb decreases with time on 
account of the bubble's inflation and uprising, the actual situation is close to the case n 2 < 1, which explains why the system 
with the bubble's deceleration included tends to live slightly less than the other. 

The solid line plots the perturbation's evolution as described by the full equation (40). The amplitude oscillates, and 
initially decays as <j> oc t~ 1//10 on account of the "Hubble drag" (discussed after Equation (47)). After several equivalence 
times t cq the Rayleigh-Taylor instability finally takes over and disrupts the bubble. Note that the effect of the Hubble drag 
is fundamental in lengthening the life time of a bubble. Incidentally, note that a very similar effect also occurs in the context 
of the evolution of cosmological perturbations; the Hubble expansion inhibits the perturbations, and in particular the density 
contrast grows as a mere power-law of time (e.g. Peacock 1999). 

Figure 3 shows the results when the the initial perturbation is cj> = 10 -2 . All other parameters as in Figure 2. 

In Figure 4 the wavelength is one third that in Figure 2; all other parameters are the same. As expected, perturbation 
growth faster. 

In Figure 5 the results of releasing the bubble at r c = 15 kpc are presented. The bubbles lives a little longer. This will 
be discussed below. 

Changing the initial condition of 4> does not change much the results, as is shown in Figure 6; all other parameters as in 
Figure 2. 

Figure 7 shows the effect of starting the perturbation at t — 10 6 yr instead of t = 10 3 yr; all other parameters as in 
Figure 2. 

Finally, in Figure 8 we plot the modulus of the relative velocity v$ = Rb<j> of the perturbation with respect to the bubble 
front; note that at early times the perturbation is oscillating: this is the over-stable regime for the three cases calculated in 
Figure 2, which includes bubble acceleration, with (solid line) or without (dot-dashed line) the drag term. When the bubble 
acceleration is considered without the drag term (dot-dashed line), the perturbations attain quite high velocities. In the inner 
regions of cooling flow clusters the sound speed is 

~keV ) ' ^ 
and therefore the perturbations oscillate with Mach numbers of A4 ~ 0.1 — 0.2. At such fast motion some dissipation, which 
is not included in our calculation, is expected. The dissipation will slow down the growth of the instability in the over-stable 
regime (the oscillatory phase), making the bubble more stable already when the drag term is not included (dot-dashed line 
in Figures 2 and 8). 

From the Figures some robust trends emerge. 

(i) Life time. The life time of the bubble is always significantly (up to an order of magnitude) longer when the initial 
inflationary phase is accounted for (Soker et al. 2002). For typical values used here, the bubbles are torn apart by the RTI 
after a time > 10 s yr, assuming that the jet is active for ~ 3 x 10 7 — 10 s yr (Birzan et al. 2004). This is enough to account 
for the existence of outer ghost bubbles in clusters of galaxies, whose ages are < 10 8 yr (Dunn et al. 2005). 

(ii) Bubble's location. The bubble injected close to the cluster centre live shorter than those injected at larger distances. 
The reason is that the growth time of the RTI is proportional to Xj/ 2 (Equation (1)). Bubbles near the centre expand at a 
lower rate (see Equation (3)), hence for a constant ratio khRb = 2nRb/\h bubbles near the centre have shorter perturbation 
wavelength Xh, and therefore the RTI evolves faster. Also, near the centre gravity is higher, also shortening the RTI growth 
time. Note that the effect of gravity is not very large in the Navarro et al. (1997) gravitational potential we have considered, 
but would be more relevant in presence of a steeper gravitational profile. 
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(iii) Wavelength. The instabilities with short wavelength (i.e. higher value of yc = k^Rt) develop first (see also Equation (1)). 
They may change somewhat the bubble boundary, but it is the large wavelength perturbations, ~ that will tear the 
bubble apart. 

(iv) Energy Injection Rate We have run some calculations with a lower energy injection rate E, without reporting the 
relative figures, since they are very similar to those already presented. As one might expect, a lower energy injection rate 
entails a weaker deceleration and a weaker drag, and hence a shorter bubble lifetime. For E = 10 42 erg s _1 , for instance (i.e. 
two orders of magnitude less than the case presented in the paper), the equivalence time only halves: t cq ~ 10 7 yr instead of 
t cq ~ 2 x 10 7 yr. The difference is not large, and is due to the weak dependence of the bubble radius on the energy injection 
rate (Rb ex. E 1 ^ 5 , see Equation (3)). 

Before ending this Section, it is worth to compare our results with the numerical simulations of Briiggen et al. (2002), 
who studied numerically the evolution of radio bubbles during the active inflation phase. A detailed comparison is difficult, 
since Briiggen et al. (2002) include the effects of a possible Kelvin-Helmholtz instability (KHI), which we have neglected. 
The KHf resulting from the bubble shear velocity with respect to the outer environment shapes the bubble as a mushroom, 
especially for low energy powering AGN jets. The RTI, on the other hand, tends to disrupt the top of the bubble first. As it is 
apparent from Figures 2, 4 and 5 of Briiggen et al. (2002), the top of the bubble keeps its round shape even if the jet is weak, 
showing that the RTI is not very relevant at this stage. Also the dependence of the lifetime of a bubble on the E supplied 
by the jet agrees qualitatively with our results, although a more quantitative assessment is difficult. We also note that the 
flow of the ICM along the bubble sides as the bubble rises, will stretch the magnetic field lines along the flow direction. This 
builds the right tension to suppress the KHI (De Young 2003). It is much more difficult to build a magnetic field to stabilise 
the top front against the RTI, however. Therefore, although we claim than magnetic fields are not required to protect the 
bubble against the RTI, we do accept that they can suppress the KHI. 

6 THE POST-INFLATION PHASE 

The results presented in the previous Sections refer to an early phase of the lifetime of a bubble, when it is still powered 
by the jet launched by the central AGN. It is important to stress that this is a crucial assumption for our model. Yet, the 
jet cannot power the expansion indefinitely, and after few 1 7 yr (e.g Alexander and Leahy 1987; Birzan et al. 2004) the jet 
fades. At the end of this "active" inflation phase the amplitude of the perturbations has been slightly suppressed by the drag 
and deceleration effects (see Equation (45), or Figure 2), so the instabilities will take more time to develop fully. After the jet 
switches off, the bubble rises buoyantly due to its low density. The process is almost adiabatic (there are but negligible heat 
exchanges between the bubble and the environment), so the bubble's expansion is far too slow to induce a deceleration or a 
drag force strong enough to hinder the growth of the RTI. The RTI sets in, leading to the eventual disruption of the bubble 
itself. 

At this stage the evolution of the bubble is quite difficult to follow analytically, on account of a host of new dynamical 
effects which will be discussed below. We refer to the existing literature for the numerical simulations of the evolution of these 
buoyant bubble, no more powered by an active jet (we refer in particular to Briiggen and Kaiser 2001 and Reynolds et al. 
2005: since they do not include magnetic fields in their simulations, their results are more directly comparable to the case 
analysed in the present work). 

In this Section we address briefly the most relevant physical effects occurring in the buoyant phase, as well as their 
influence on the growth of the bubbles' instabilities. 

After the end of the inflation's stabilising effect (or after few times t e q, when the deceleration and the drag are too weak 
to be relevant) the RTI is free to evolve. As the perturbations' amplitude attains the magnitude </> p h ys /Ah ~ 0.1 — 0.2, i.e. 
<j) ~ 0.5 — 1/x the RTI enters the non-linear regime. Ofer et al. (1992) have shown that if the RT perturbations are not exactly 
monochromatic (which is the case in a realistic situation) the coupling between modes of different wavelengths suppresses the 
fast-growing short-wavelength modes; the overall effect is to suppress the growth of the instabilities, resulting in a further 
delay in the bubble disruption. 

During the non-linear evolution the relative velocity of the bubble with respect to the ambient ICM may significantly 
affect the dynamics, with potentially dominant contributions from the KHI. There is another reason why the KHI becomes 
important only at this stage: in the early inflation phase the density contrast between the bubble interior and the ambient gas 
is very large, and in this situation the KHI is suppressed (Kaiser et al. 2005). Later on, however, when the density contrast 
lowers, the KHI may become important enough to be the primary cause for the bubble disruption. 

Some authors suggest that the KHI growth rate may be reduced or suppressed in this phase. Kaiser et al. (2005) and 
Reynolds et al. (2005) show that the transport phenomena (plasma viscosity and thermal conduction) are able to suppress 
the growth of the fastest Rayleigh-Taylor and Kelvin-Helmholtz modes. This result is particularly important in the context 
of the present paper, since in the absence of an ordered magnetic field parallel to the bubble surface neither viscosity nor 
thermal conduction are expected to be strongly suppressed. 
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In the different context of the stability of cold fronts against the KHI Churazov and Inogamov (2004) have shown that 
if the interface between the interior of the bubble and the ICM has a finite thickness, the fast short-wavelength KHI modes 
are suppressed. As mentioned before, we do not expect the magnetic field to suppress the RTI, but we do expect it is able to 
reduce the growth of the KHI modes. 

To summarise, even in the buoyant phase the growth of the instabilities may be reduced by the occurrence of several 
effects. We repeat that we do not expect them to prevent the final disruption of the bubble. The bubble will be torn apart, 
but on time scales of few 10 8 yr, consistent with the estimated ages of the "ghost" bubbles observed in several cooling flow 
clusters. The final onset of the RTI, with consequent disruption of the bubble, is only delayed but not avoided, which is 
consistent with ghost bubbles showing RTI features. Soker et al. (2002) suggest that the protrusion in the northwest bubble 
of the Perseus cluster (Fabian et al. 2000, 2002) is a clear signature of the late onset of the Rayleigh-Taylor instability. 

7 SUMMARY 

To summarise, our main conclusion is that the life time for X-ray deficient bubbles in cluster of galaxies, and the development 
of RTI features in ghost bubbles can be explained from pure gas dynamical effects, and there is no need to assume the 
existence of stabilising magnetic field. This conclusion is robust, and is not affected by the several simplifying assumptions we 
have made. For typical parameters, inclusion of the bubble inflation phase increases the expected life of bubbles by a factor 
of 5 - 10. 
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Figure 1. Upper panel: evolution of the radius of a bubble R b (solid line) compared to the distance r c between the centre of the bubble 
and the centre of the cluster (dashed line). Middle panel: expansion velocity R b of the bubble front. Lower panel: acceleration R b of the 
bubble front. The dotted horizontal line shows the modulus \g\ of the gravitational acceleration at the position where it equals \R b \ on 
the top of the bubble. In this example the bubble is released at the initial distance r c = 5 kpc from the centre. 
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Figure 2. This plot compares the evolution the amplitude of a RTI of a perturbation in three cases. Plotted is the absolute value of the 
co— moving perturbation <f>, i.e. the ratio between the real physical perturbation amplitude and the bubble's radius Ri, (Equation (41)). 
The cusps are artifacts of the logarithmic scale. In this example the bubble is released at r c = 5 kpc from the centre, and the ratio 
between the circumference of the bubble and the wavelength of the physical perturbation is >c = 8. The perturbation starts 10 3 yr after 
the bubble is born, with an initial co-moving amplitude </>q = 10~ 3 and zero relative velocity <j>o = 0. The dashed line plots the case in 
which the inflation's deceleration is neglected (Equation (58)). The dot-dashed line plots the evolution of the perturbation without the 
"Hubble drag" (Equation (59)). The solid line is the evolution described by the full equation (40). The dotted vertical line marks the 
time t cq ~ 2 X 10 7 yr when the acceleration R\ of the bubble front equals in magnitude the gravitational acceleration on the top of the 
bubble: \Rb\ = \g(Tc + Rc)\- After 10 s yr the RT perturbation is well developed in the cases described by Equations (58) and (59), but 
it is still small in the case described by Equation (40). 
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Figure 4. The lines have the same meanings as in Figure 2; also the parameters arc the same, except the factor x = k^Rf, = 24. 
The short wavelength perturbations grow faster (see Equation (1)), although they are not responsible for the final RT disruption of the 
bubble. 
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Figure 6. In this graph we show the effect of different initial velocities on the perturbation's amplitude evolution, as described by 
Equation (40). The bubble is released at r c = 5 kpc from the centre. The integration starts 10 3 yr after the inflation starts, with a 
parameter x = 8 and an initial relative amplitude <f>0 = 10 -3 . The dashed line corresponds to the initial condition (57b) with t p = t b /l0, 
the solid line to t p = t^, and the dotted-dashed line to t p = 10 tf,. As expected, the initially fastest-growing perturbations develop 
first, although the effect is not very large. The perturbations reach the amplitude \<f>\ = 0.5 (at which we may considered the bubble as 
disrupted by the instability) at the times t ~ 3.3 X 10 8 yr (case t v = i;,/10), t ~ 6.8 X 10 s yr (case tp = %), t ~ 1.8 X 10 9 yr (case 
tp = 10 tt,); there is a factor ~ 6 difference, while the ratio ip/tf, spans two orders of magnitude. In all cases, the perturbation develops 
less than in the case without inflation (compare with Figure 2). 
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Figure 7. In this graph we show the effect of the choice of the initial integration time on the perturbation's amplitude evolution, 
as described by Equation (40). The parameters and the meaning of the lines are the same as in Figure 2, but the integration starts 
to = 10 6 yr after the inflation onset. The difference between the case without acceleration (dashed line, Equation (58)), and the case 
without the Hubble drag (dot-dashed line, Equation (59)) is smaller. In the case without acceleration the repulsive force x\g\/ R^t) is 
smaller, since Ri, starts from a larger value. In the case without the Hubble drag the oscillations do not have enough time to develop 
large enough amplitude before the time t eq when they become unstable. 
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Figure 8. In this graph we show the modulus of the relative velocity = cj> of the perturbation with respect to the front of the 
expanding bubble. The parameters are the same as in Figure 2. 



